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ABSTRACT 

When  fine  wires  are  exploded  by  discharging  a  capacitor, 
the  circuit,  exhibits  damped  electrical  oscillations  as  an  U-R-C 
circuit;  however,  the  resistance  varies  during  the  explosion 
which  complicates  the  analysis*  To  date  the  resistance  has  been 
calculated  either  by  integrating  twice  in  succession  the  record 
of  the  rate  of  change  of  the  current* ,  or  by  the  current- voltage 
method^.  A  different  approach  is  proposed  in  which  the  curcuit 
equation  is  differentiated  with  respect  to  time  and  then  solved 
by  a  method  similar  to  the  WKB  method.  The  method  will  be 
discussed,  original  data  provided,  and  values  of  calculated 
resistances  obtained  by  the  double  integration  method  will  be 
used  to  indicate  the  validity  of  the  solution. 
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I.  introduction 

A*  evidenced  by  publication*  in  the  literature,  conferences  *uch  as  the 
present  one,  and  excellent  adaptation*  to  industrial  fields,  experimenter*  have 
been  successful  in  using  exploding  wires  for  many  different  task*.  They  have 
done  thi*  with  a  minimum  of  theoretical  work  and  a  maximum  of  experimental 
effort*  Thi*  work  represent*  another  attempt  to  combine  empirical  measure¬ 
ment*  with  basic  theory  to  determine  whether  a  postulated  resistance  varia¬ 
tion  can  fit  the  observed  data. 

We  shall  solve  the  circuit  equation  with  a  constant  resistance  as  a  model; 
solve  the  same  equation  with  a  resistance  that  varies  with  time  using  a  method 
similar  to  that  used  by  Wenzel,  Kamera,  II  Brillouin  to  indicate  the  algebraic 
form  required  for  the  resistance;  use  this  form  to  establish  the  terms  of  a  ser¬ 
ies  solution  for  the  current  in  the  circuit;  and  lastly  show  how  the  solution  can 
be  made  to  fit  the  observed  data. 


II.  MATHEMATICAL  PROCEDURE 

The  basic  circuit  equation  for  the  capacitor  discharge  circuit  is  (as  usual) 


t 

l-5T  +ri  *  ir  S  m  -  vo 


wh^re  L  is  the  inductance  of  the  circuit,  I  is  the  current  in  the  circuit,  t  is  the 
time  measured  from  switch  closure,  R  is  the  sum  of  the  wire  resistance  and 
the  resistance  of  the  circuit  element*,  C  is  the  capacitance  of  the  circuit,  and 
V  is  the  voltage  to  which  the  capacitor  was  charged  at  switch  closure  time* 
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To  solve  this  equation.  It  Is  first  differentiated  with  respect  to  time  to  eliminate 

the  integral.  In  this  wor!<  the  L,  C,  and  VG  are  assumed  to  be  independent  of 

4  * 

time.  The  equation  becomes 


L 


i  R 


dl 

dt 


+  I 


dR 

dt 


(2) 


A.  Solution  with  R  Constant 

At  the  risk  of  repetition  we  will  solve  this  equation  with  a  constant 
resistance  to  form  a  limiting  case  for  a  later  solution.  One  assumes  a  solu¬ 
tion  with  two  constants  A  and  B  which  are  evaluated  with  the  initial  conditions. 


* 

T  .  wt  ^  w  t 

1  »  Ae  +  Be 


ay 


w  and  w5,i  are  taken  to  be  complex  conjugates  with  A  the  real  part  and  B  the 
imaginary  part. 

\v  *  a  +  ib  (4) 


The  assumed  solution  is  set  into  the  differential  equation.  By  setting  the 
coefficients 

Ao^tLw2  +  Rw  +  -L)  +  De'V  *  {L*2  +  R%v*  +  -L  )  .  0 


(5) 


•£  the  exponentials  independently  equal  to  £erO  one  solves  to r  a  ana  b  In 
terms  of  L.,  H#  and  C*  This  leads  to 


a  = 


R 

2L 


16) 
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Thus, 
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With  the  usual  initial  conditions  that 


A 


I 


Idt  =  0 


t  =  0 


(9a) 

(9b) 


the  current  is  fovnd  to  be 


(10) 
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B*  Solution  with  R  a  Variable  in  Time 


The  WKR  method  of  solving  this  type  of  equation  assumes  that  the 
variation  of  R  is  small  with  respect  to  itself  and  takes  place  in  a  time  short 
compared  with  the  basic  period  of  oscillation.  While  It  does  vary  during  a 
small  portion  of  the  base  wave  length  in  the  solution,  it  varies  by  a  large 
amount;  therefore,  we  shall  not  be  able  to  follow  the  WKb  method  completely 
to  simplify  the  differential  equation.  However,  We  shall  use  the  general  pro¬ 
cess  to  obtain  an  equation  for  solution. 

As  before,  we  assume  a  solution  Wt  allow  the  w#S  (which  means  the 
a*s  and  the  b’s  also)  to  be  a  function  of  time.  In  this  cafe,  with  **bitraty 
constants  C  and  D  intended  to  be  determined  by  the  initial  conditions,  we  have 

(it) 


in  which  the  variation  ofw  with  time  is  inserted  by  employing  the  integral 
representation.  This  solutionis  substituted  Into  the  euqatioi>  to  obtain 
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As  before,  we  substitute  in  the  original  equation  with  w 


a  +  ib;  set  the 


coefficients  of  the  exponential  terms  separately  equal  to  zero;  separate  the 
real  and  imaginary  parts  of  these  equations  and  solve  for  a  and  b  in  terms 
of  R,  E.  C  and  t.  This  process  leads  to 


1  jib^  R 
1  "  ‘  "2b  "dt-  '  2L 

2  2—2 
1  d  b  3  /  db  \  2  _  dt  R  1 

2b  2  ”  2  \  dt  /  2L  '  _  2  LC 

dt  4b  4L 


(13) 


(14) 


The  equation  for  b  obtained  by  using  w*  =  a  -  ib  will  ba  th(j  same  as  Equation 
14. 

At  this  point  in  the  WKB  process,  the  second  derivatives  and  products 
of  first  derivatives  are  assumed  to  be  small  enough  to  be  neglected.  In 
Equation  14  the  solution  for  b  would  be  at  hand;  hoWeVar,  given  the  solution 
for  b  in  that  case  the  term  with  the  first  derivative  Is  nO  longer  neglectable. 

If  one  includes  this  term  and  solves  for  b,  he  finds  that  the  value  of  the  second 
derivative  term  is  now  no  longer  neglectable.  Therefore*  no  terms  are 
neglectable  in  this  case.  The  complementary  solution  of  I^qtiation  14  Is  found 
to  be 

•  • 

2U+7+T7)  ‘  2(t+r  -  is)  -  (t+  r)2  ~  (1  ' 
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in  which  r  and  s  are  the  two  undetermined  constants  required  for  a  second 
order  differential  equation. 

The  particular  solution  of  this  equation  depends  upon  the  form  of  the 
right  hand  side  since  it  is  the  driving  function.  Usually  a  solution  is  rela¬ 
tively  easy  to  obtain  if  the  driving  function  is  of  the  same  form  as  the  comple¬ 
mentary  solution.  Therefore,  we  let  R  be  expressed  in  the  form  of  a  truncated 
series.  Then,  b  may  be  obtained  as  a  series  of  this  same  form  that  might 
hopefully  be  truncated  without  a  great  loss  of  accuracy.  A  series  for  R  with 
terms  similar  to  Equation  15  if  s  is  suppressed  in  favor  of  r  might  be 


R  =  k  + 


(t  +  g) 


(t  +  g) 


(t  +  g) 


as  an  equation  in  which  k,  h,  gf  f,  e  and  g  are  constants.  We  assume  accord- 
ingly  a  series  solution  for  b 


s  m  . 

b  =  2  - (17) 

j=°  (t  +  n) J 


set  it  into  the  differential  equation,  set  the  coefficients  of  (t  +  n)  '  terms 
separately  equal  to  zero  and  determine  the  coefficients  mj. 

After  a  series  of  algebraic  manipulations,  b  is  found  to  be 


(IB) 
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and  a  is  determined  as 


k 

h 
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in  which 


W  = 


1 

LC 


41, 
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S  =  -  2 
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(24) 


The  first  three  terms  of  the  exponent  ( 


kt  h  t  +  g 

E  =  '  TT'TT  ln - ~ 

r  2L  2L  g 


wdt)  are  given  as  a  real  part 


“/  h  \  /  U  \ 

/■ 

+  -i- 

t 

2W2 

+  2  L 

g(t  +  g) 

(25) 


and  an  imaginary  part 
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2  W 


2  W 


+  g) 


(26) 


Given  the  initial  conditions  and  the  terms  Er  and  E^,  one  finds  that  the  expres- 
sion  for  the  current  is  as  could  be  shown  in  part  A 


C.  Comparison  of  Solutions 

The  following  statements  may  be  made  about  the  solution  for  the 
constant  and  variable  cases: 

(1)  The  solution  for  the  variable  R  case  being  a  series  should  reduce 

to  the  constant  R  case  in  the  limit.  We  find  this  to  be  so:  when  h  =  g=  f=  e  =  0, 
Er  in  Equation  25  reduces  to  Equation  6  and  E^  reduces  to  Equation  7. 

(2)  From  experiments,  the  current  approaches  a  damped  oscillation 
which  implies  fairly  constant  resistance.  We  find  that  R  approaches  a  steady 
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value  given  by  k  as  the  time  gets  large.  Looking  at  Er,  we  find  that  the  third 

(and  succeeding  although  not  shown)  terms  approach  a  steady  value  because 

— —  approaches  one.  The  second  term  continues  to  increase  but  not  as 
t  +  g 

rapidly  as  the  first  term.  Thus,  the  curve  decays  at  a  slightly  larger  rate 
for  the  variable  than  for  the  constant  resistance. 

exhibits  the  same  behavior.  However,  since  E^  is  the  argument  of 
the  sine  function,  the  main  effect  is  to  produce  a  modulation  of  the  phase  or 
frequency.  Thus,  one  wd  uld  expect  a  large  variation  at  the  beginning  of  the 
discharge  and  slower  one  at  the  end  of  the  discharge.  These  changes  will 

i 

cause  the  curve  to  be  nonsinusoidal. 

The  amplitude  factor  contains  dE^/dt  which  is  also  a  function  of  time. 
Initially  the  variation  is  large  but  will  damp  out  because  the  time  varying 
terms  have  t  in  their  denominator.  Thus,  for  long  times  we  find  that 

dE. 

1 

~dt" 

which  is  found  in  the  solution  for  a  constant  R.  (See  Equation  10.  ) 

III.  EXPERIMENTAL  DATA 

The  above  expression  for  the  current  should  be  fitted  to  the  data  so 
that  the  various  constants  might  be  obtained.  Because  the  various  terms 
are  made  up  of  a  series  expansion  such  curve  fitting  would  be  quite  compli¬ 
cated.  In  place  of  this  method  we  shall  use  the  double  integration  method  to 
determine  the  value  of  R. 
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From  the  original  circuit  equation  the  resistance  is  given  as  shown 


in  reference  1. 


(29) 


Since  the  rate  of  change  of  the  current  is  measured  with  a  Rogowsky  loop,  the 
current  is  obtained  by  integrating  once  and  the  integral  of  the  current  by  inte¬ 
grating  twice.  With  the  calibration  constant  one  can  obtain  the  value  of  the 
resistance  as  shown  in  the  equation.  If  one  uses  a  current  shunt,  the  other 
terms  are  obtained  by  differentiating  and  integrating  the  data  for  the  current. 

Since  every  experimenter  has  different  equipment  and  since  his  results 
are  largely  a  matter  of  his  instrumentation,  Figure  1  shows  the  exploding  wire 
apparatus  used  to  obtain  our  data.  The  12  capacitors  can  be  charged  to  50  KV 
thereby  having  a  storage  capability  of  10^  joules.  The  ringing  frequency  is 
275  kilocycles  so  that  the  pulses  are  made  up  of  3  or  4  cycles  that  carry  the 
bulk  of  the  power.  The  wire  commonly  used  is  2  mil  tungsten  1/2  inch  long. 
Figure  2  shows  the  cartridge  chamber  housing  the  wire. 

A  Rogowsky  loop  measures  the  rate  of  change  of  the  current  as  shown 
in  Figure  3.  The  sweep  speed  of  the  Tektronix  scope  is  2  microseconds/cm 
which  speed  is  used  to  obtain  the  values  of  L  and  R  in  the  main  part  of  the 
pulse.  Figure  4  shows  the  pulse  taken  at  1/2  microsecond/cm.  This  shows 
the  shape  of  the  first  half  cycle.  After  reading  the  deflections  with  a  12X 
microscope,  data  from  the  curve  is  plotted  as  shown  in  Figure  5. 
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The  frequency  is  determined  from  the  cross  over  points  thereby  giving 


p 


the  inductance  from  the  formula 


si|j _ si 

2ff  V  LC  .  -  2 
4  L 


(30) 


Using  curves  as  shown  in  Figure  3  the  value  of  R  in  the  steady  portion  of  the 
decay  curve  is  found  to  be  about  0.030  ohms.  Using  this  value  of  resistance 
the  inductance  is  found  to  0.111  x  10”  ^  Henries.  Using  these  numbers,  numerical 
integration  gives  us 


I 


and 


(31) 


(32) 


from  which  the  circuit  resistance  is  obtained  for  each  interval  of  time.  This 
relation  is  also  plotted  on  Figure  5. 

By  standard  curve  fitting  methods  the  equation  for  the  resistance  curve 
may  be  obtained,  when  t  is  expressed  in  microseconds  and  R  in  ohms. 


R 


0.  030 


0.  0244 

t  +  o.  61 


0.  048 

(t  +  0.  oi)2 


0.  00650  0.  000271 

(t  +  0.  01)3  (t  +  0. 01)4 


(33) 


This  shows  that  the  resistance  is  approaching  the  limit  0.  030  as  time  increases. 
It  also  indicates  that  many  terms  are  needed  to  get  an  accurate  value  because 


1  1 


the  coefficients  do  not  decrease  rapidly.  This  indicates  that  many  terms  will 
be  needed  in  the  expression  for  the  current. 

To  obtain,  a  possible  rationale  for  the  particular  form  of  the  equation 
for  the  resistance ,  one  might  propose  the  following  derivation.  According  to 
Lin°  (and  Bennett')  the  radius  of  the  shock  wave  proceeding  outward  from  a 
line  explosion  is  proportional  to  the  square  root  of  the  time  since  the  explosion. 
If  one  thinks  of  the  current  path  as  being  the  material  behind  or  within  that 
shock  wave,  then  the  cross  sectional  area  will  be  proportional  to  that  radius 
squared.  By  the  simplest  theory  of  conductors  the  resistance  is  inversely 
proportional  to  the  cross-sectional  area.  Therefore,  the  resistance  of  the 
circuit  will  have  a  constant  term  for  the  busbars  and  other  items,  and  it  will 
have  a  term  for  the  wire  or  plasma  or  gap  that  is  inversely  proportional  to  a 
time  interval.  Since  the  resistance  is  also  a  function  of  the  number  of  current 
carriers,  their  mean-free -path,  and  their  temperature  or  internal  energy,  it 
is  not  surprising  that  other  terms  should  be  included  in  the  series  expression 
for  the  total  resistance. 


IV.  CONCLUSIONS 

(1)  The  circuit  equation  was  solved  by  the  series  method.  (2)  A  pro¬ 
cedure  was  shown  that  indicated  the  type  of  series  that  should  be  used.  (3) 
This  series  was  found  to  be  amenable  to  the  data,  but  a  series  was  obtained 
that  converged  rather  slowly.  (4)  Terms  were  found  in  the  series  which  would 
represent  variations  in  frequency  and/or  phase. 
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When  fine  wires  are  exploded  by  discharging  a  capa¬ 
citor,  the  circuit  exhibits  damped  electrical  oscillations 
as  an  L-R-C  circuit;  however,  the  resistance  varies  during 
the  explosion  which  complicates  the  analysis.  To  date  the 
resistance  has  been  calculated  either  by  integrating  twice 
in  succession  the  record  of  the  rate  of  change  of  the  cur¬ 
rent,  or  by  the  cur  rent- voltage  method.  A  different 
approach  is  proposed  in  which  the  circuit  equation  is  ciif  — 
derentiated  with  respect  to  time  and  then  solved  by  a  me 
method  similar  to  the  WKB  method.  The  method  is  dis¬ 
cussed,  original  data  provided,  and  values  of  calculated 
■resistances  obtained  by  the  double  integration  method  are 
used  to  indicate  the  validity  of  the  solution. 
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